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Abstract 

Evolution PDEs for dispersive waves are considered in both linear and 
nonlinear integrable cases, and initial-boundary value problems associated 
with them are formulated in spectral space. A method of solution is pre- 
sented, which is based on the elimination of the unknown boundary values 
by proper restrictions of the functional space and of the spectral variable 
complex domain. Illustrative examples include the linear Schrodinger 
equation on compact and semicompact n- dimensional domains and the 
nonlinear Schrodinger equation on the semiline. 



1 Introduction 

Initial-Boundary Value (IBV) problems for Partial Differential Equations (PDEs) 
play an important role in applications to physics and, in general, to natural sci- 
ences. 

It is well known that the basic difficulty associated with the study of IBV 
problems for linear and nonlinear integrable PDEs is the presence of unknown 
boundary values in the relevant equations of any method of solution. In this 
paper, after formulating the IBV problems for linear and nonlinear integrable 
PDEs in spectral space, we present a method of solution, the Elimination-by- 
Restriction (EbR) approach, which is based on a strategy of elimination of the 
unknown boundary values by proper restriction of the functional space and of 
the complex domain of definition of the involved spectral functions. This ap- 
proach is inspired by the Green's function (GF) method, which, in the linear 
context, is essentially its counterpart in configuration space, and by our recent 
findings in the nonlinear context. The paper is organized as follows. In §2 we 
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deal with IBV problems for linear PDEs with constant coefficients. After in- 
troducing the proper Fourier Transform (FT) for that problem and establishing 
its analyticity properties, we first express the Fourier transform of the solution 
in terms of the Fourier transforms of known and unknown initial - boundary 
values using Green's formula. Then we present the EbR approach in which, 
using systematically a strategy of elimination of the unknown boundary values, 
one obtains the appropriate spectral representation of the solution whose sup- 
port may eventually turn out to be discrete rather than continuous as in the 
general Fourier integral one starts from. We illustrate the power of the method 
solving IBV problems for the Schrodingcr equation on an n-dimcnsional rectan- 
gular box and quadrant. In §3 we apply the EbR strategy to soliton equations. 
After defining the proper spectral transform S(k,t) for the given IBV problem, 
we apply the EbR procedure eliminating the unknown boundary values from 
the equations defining S(k, t) and characterizing S(k, t) via a nonlinear integral 
equation. This approach, presented on the prototype example of the nonlinear 
Schrodingcr (NLS) equation on the semifine, can be in principle generalized to 
the segment case. Since the EbR technique operates on the spectral variables 
conjugated to the space ones, it works well either using the space-time trans- 
form or just the space transform. In the linear case we find it simpler to work 
with the space-time FT, while in the nonlinear case it seems more convenient to 
use the space transform, namely the well-known Inverse Scattering (Spectral) 
Transform (1ST). We finally show the equivalence between IBV problems for 
soliton equations on the semifine and some specific forced initial value problems 
on the whole line. An important application of this equivalence is that, from the 
well-known asymptotics of soliton equations on the whole line, one can obtain 
immediately the asymptotic behaviour for IBV problems on the semiline with 
decaying boundary data. We have confined the relevant literature to §4. 

The results contained in this paper are an expansion of part of the material 
presented by the authors at the Euroconference "NEEDS 2001" and at the 
workshop: "Boundary Value Problems" , the first and last events of the Semester: 
Integrable Systems, held at the Isaac Newton Institute of Cambridge during the 
period July - December 2001. 

2 The Elimination-by-Restriction Approach: 
The Linear Case 

It is well -known that the Fourier Transform (FT) is the proper tool to solve 
initial - boundary value (IBV) problems for linear PDE's in lZ n+1 with decaying 
boundary values: 



£(V,|Wx,i) = /(x,t), 
u(x, 0) = u (x), 



x = (xi, ..,x n ) G TZ n , t > 0, 
u(x, t) — ► 0, |x| — ► oo, 



(1) 
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where y = ••, gf - ); £ is a constant coefficients partiai differentiaf operator, 
u(x, £) is the unknown field, /(x, i) is a given forcing and «o(x) is the given initial 
condition. 

In this section we present a very effective approach, in Fourier space, for 
solving more complicated IBV problems, defined on compact or semi - compact 
space domains V: 

£(V,^Wx,t)=/(x,t), xeRK", <>0, (2) 
with Dirichelet or Neumann or mixed boundary conditions on dV . 

2.1 The Fourier Transform and its properties 

The natural FT associated with the space - time domain V — V ® (0, oo) (in 
short: FT-p) is defined by 

F(k,g)= / rfxdte- l(k ' x+9t) F(x,i) (3) 
Jv 

for any smooth function F(x, i), (x, t) 6 I?, assuming that F(x, t) — * 0, t — > oo 
fast enough; here k = (fci, .., k n ) £ 7?. n , q & TZ and k • x = fcjXj. Its inverse: 

F(x,t) X p(x,t) = ^ ^.e^-^FOc,?) (4) 

reconstructs F(x.,t) in D and zero outside, where Xx>(x, t) is the characteristic 
function of the domain V: Xz>(x, t) = 1, (x, t) G I?, xr>( x , t) = 0, (x, t) £ V 
(therefore: XT>(*,t) = Xv(x-)H(t), where H(t) is the usual Heaviside (step) 
function). 

If the space domain is the whole space: V = lZ n , the FTt> (§) is defined in 
A = 1Z n ® T q , where 2 q is the closure of the lower half q-plane 2 q , analytic in 
q &I q , Vk 6 lZ n and exhibits a proper asymptotic behaviour for large \q\ in the 
analyticity region. If the space domain V is compact, the FTx> acquires strong 
analyticity properties in all the Fourier variables: it is defined in A = C n <S> T q , 
analytic in q £ I,, Vk £ C n , entire in every complex kj, j — 1, ..,n V5 6 I, 
and exhibits a proper asymptotic behaviour, for large (k, q), in the analyticity 
regions. If the space domain is semi - compact, then the analyticity in the 
Fourier variables kj, j = 1, .., n is limited to open regions of the complex plane, 
depending on the geometric properties of the domain V. 

To express the FTx> of the solution in terms of the FTjy's of the forcing 
and of the IB conditions we make use of the well - known Green's formula 
(identity): 

bCa — aCb = div J(x, t), (5) 
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and of its integral consequence, the celebrated Green's integral identity: 

/ (bCa — aCb)dxdt = / J(x, t) ■ vda, (6) 
Jv JdV 

obtained by integrating over the domain T> and by using the divergence 
theorem. In equation (||), C is the formal adjoint of L: £ = C(—S/, — Jj), >^( x , t) 
is an (n + l)-dimensional vector field, div is the (n + l)-dimensional divergence 
operator and a(x, t) and &(x, t) are arbitrary functions. In equation (||), da is 
the hypersurface element of the boundary and v is its outward unit normal. We 
remark that, given C, its formal adjoint C and two arbitrary functions a and 
6, an (n + l)-dimensional vector field J(x, t) satisfying the Green's formula (|J) 
always exists and can be algorithmically found to be a linear expression of a, b 
and their partial derivatives of order up to N — 1, if C is of order N . 

The arbitrariness of a and b allows one to extract from (|^) and (|^) several 
important informations on the BV problem; with the particular choice 

a = u{x,t), b = e-^ x+ rt/£(ik,iq), (7) 

where £(ik, iq) is the eigenvalue of the operator C, corresponding to the eigen- 
function e^ kx+<3t - ) , the vector field J takes the following form: 
J = e ~ , ( k ' x +9 t ) J'(x, t; k, q)/C{ik, iq) and the Green's integral identity (||) gives 
the FT-p of the solution in terms of the FTp's (or, maybe, of generalized FT's) 
of the forcing and of all the initial - boundary values: 

ft(k; q) />. g) - lev 6 - (k - X+?t) ^(x, t; k, g) • uda _ : ^(k, g ) 

£(ik, zg) £(ik, ig)' ' 

(8) 

In general, £(ik, ig), the denominator of the above equation, is an entire and, 
most frequently, polynomial function of all its complex variables and its zeroes 
may lie on the real axis; therefore, before calculating the inverse FT, we must 
regularize it: 

C(ik, iq) -» Creg («k, iq) ; (9) 

i.e., we must move a bit the singularities off the real axis, outside the domain 
A. 

Its inverse transform (|J) gives the corresponding Fourier representation 
of the solution: 

U M - »(x, t )x„(x, t ) . f ^e'<— (x, t ) e K-» 

(10) 

Clearly this is not the end of the story since, in general, the RHS of equation 
(§h depends on known and unknown boundary values. 
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2.2 Elimination-by-Restriction in Fourier space 

The traditional ways in which IBV problems for linear PDE's are solved consist 
in finding convenient strategies for eliminating the unknown boundary conditions 
from the representation of the solution. On this idea is based the celebrated 
Green's function (GF) approach, in which: 
i) one constructs the Green's integral representation 



of the solution of the IBV problem (g) as another application of (|6|) , correspond- 
ing to the choice: a(x, t) = u(x, t) and 6(x, t) = <7(x', t'; x, t), where g is any 
Green's function of L x : £ x <?(x',t';x, t) = S(x — x')8(t — t'), (x, t), (x',i') 6 D. 
ii) One uses the arbitrariness of g and constructs that particular Green's func- 
tion which allows one to eliminate contributions depending on unknown bound- 
ary values. On the elimination idea is also based the eigenfunction expansion 
method, essentially equivalent to the GF approach, in which one constructs a 
set of eigenfunctions of C with proper boundary conditions which allow again 
to eliminate the unknown boundary values. Both approaches are of functional 
analytical nature. 

In this section we shall show how the elimination strategy can be conve- 
niently implemented working in the Fourier space defined by (|^) (the elimination 
strategy in spectral space has been already investigated for soliton equations on 
the semiline (see §4)). 

Equation (Q), defined in a proper domain A C C™ +1 , usually exhibits several 
symmetry properties which are consequence of the structure of C and of the 
geometry of the space domain V. The first and more critical part of the method 
consists in constructing a linear operator £ which, by exploiting systematically 
these symmetry properties in Fourier space, annihilates, in the RHS of (p|), the 
contributions coming from the unknown boundary values: 

N 

£u(\s.,q) = £( — )(k,q) = {only known quantities}, (k, q) G A C A. (12) 

The elimination procedure, described by the linear operator £ , allows one to 
construct £u in a subspace A' of the original space of definition of u. It is 
necessary therefore to establish if this information is sufficient to reconstruct 
the solution u; i.e., if £u defines an invertible spectral transform in T>. 
Using equation (|^), this new spectral transform is defined by: 



and is invertible provided an eigenfunction ip^ q (x., t) of C exists with the prop- 




(x,<) G V. 
(11) 




(k, q) G A' 



(13) 
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erty of satisfying the completeness condition: 

^k, g (x,i)^ M (x',t')=<5(t-i')5(x-x'), (x,f) ( (x',t')€D. (14) 

In this case the inverse transform allows one to construct a function U (x, t) , 
defined in general on the whole space-time, which coincides with the solution in 
V: 

l>(x, t) = £ ^ k , g (x, t)£(£)(k, q), (x, i) e ft" ® (0, co), 

k,? (15) 

£7(x,t) = u(x,t), (x,t)e£>. 

In equations (|l4|) and ( p"5| ) £ k indicates a sum and/or an integral, depending 
on the nature of the domain A! of definition of the new spectral transform. 
Summarizing: 

Although the direct transform we started with is the Fourier trans- 
form u defined in (0), the elimination- by-restriction procedure per- 
formed on it leads to a new direct transform £u (|l3|), whose inversion 
usually differs considerably from (Q). This implies that the recon- 
structed function t/(x, t) coincides with the solution w(x, t) in V but is 
not zero outside V, inheriting the symmetry properties of the eigen- 
function <p^ q (yi, t). 



2.3 Illustrative Example 

In this section we apply the EbR approach to the IBV problem (||) corresponding 
to 



<9iA A_^n_V^9 



£ = *£ + A, A = VV=Es 

3=1 

y = {x: 0<^ <L is j = l,..,n}; 



i.e., we solve IBV problems for the (n + l)-dimensional Schrodinger equation in 
an n-dimcnsional rectangular box. 
Then: 

^ = - *J? + A ' J = {iab,b\/ a - a\/ b), , , 

£(ik,iq) = -(q + k 2 ) C reg (ik, iq) = -(q + k 2 - iO), [ 1 



where k 2 = k • k. Equations (||) and (|lj) give the following expression of the 
Fourier transform of the solution in terms of the Fourier transforms of the forcing 
and of all the initial - boundary values: 

n 

Af(k, q) := /(k, q) + m (k) + £ {[woi( k j: «) + ikjvojik^q)]- (18) 

i=i 

e -^L 3 [^.(kj, q) + ikjVLjikj^)}}, 
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in the definition domain A = C" ®I qi where /, Uo, voj, VLj, u>oj, WLj are the 
FTs of the forcing and of the following IB values: 

u (x) = u(x,t)\ t=0 , v 0j (xj,t) = u(x,t)\ Xj=0 , v Lj (xj,t) = u(x,t)\ Xj=Lj , 

^(x i ,t) = ^(x,t)| Xi=0 , WLj (^, t ) = ^,t)\ Xj=Lj ; ^ 

i.e., for instance: 

oo 

uo(k) = / dxe-* k - x Uo (x), vajO^,q) = J dt J d^e^^+^v^^t). (20) 

V Vj 

In equations (|T^)-(|2C|), as well as in the following, Xj = (x\, .., x n ) G IZ"^ 1 , 

kj = (fci, .., fcj, .., k n ) e 1Z n ~ 1 , j v dxj — J Q Ll cfcci • -(/q Lj dxj) ■ ■ J^™ dx n with the 
understanding that the superscript " indicates that the quantity underneath is 
removed. 

Of course, if all the above boundary values were known, the solution u would 
be given by the formula (|Io|): 

C/(x,<) = u(x,t) X v(x)ff(t) = -k^ ^re i(k - X+9t) ^M = 
" T0^e^+rt _ifeiL + m kn ^ ei(k .x-^^ o(k)+ 

J = lR"- 1 

e t(*jki-%|-fc a t) [i& £j -(k,-, -A; 2 ) + isignixj ~ L^Vl^, -k 2 )]}, 



(21) 



where <ikj = dk\..dkj..dk n and 7 = (ioo, 0) U (0, 00). 

In view of the distinguished parity properties of the Fourier transforms in 
in the following we shall make use of the parity operators: 

n n 

A±=n(l±o l ), A { £ = JJ(1±cti), (22) 

Mi 

where ex., is the involution <jj : fcj — > — fcj . 

Suppose we are interested in solving the Dirichelet problem; applying the 
parity operator A_ to (|q) eliminates all wq's: 



2« E [sin(%L 3 )A 0) w L j(k,-,<7) - fcj cos(fc j L J )A (j) w Lj -(k J , q)] ) /(q + k 2 - id), (k,g) £ A 

To eliminate also the wl's, the values of fcj must be restricted to the discrete 
set kj = hj :— ^ 2 -, rrij £ Z, so that the original domain A is finally restricted 



A_u(k,g) = - A_[/(k )t? ) + m (k)] +2i £ k j A ( l>v 0j (k j ,q)+ 
\ 3=1 
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to: 

7T77J 



(k, q) g A' = {(h, g); q e T q , h = (hx, ..,h„), hj = — — , nij G Z, j = 1, .., n}. 

(24) 

Therefore the EbR operator £ of this example reads: 

£• = / dk<S(k-h)A_- (25) 

and its application to u leads to the wanted result: 

£u(k,q) = (A_u)(h,g) = -{A_[/(h, g) + i«o(h)]+ 
2* £ ^AW[^(h,,g) - (-)™^ ii (h J -,g)]}/(g + / l 2 -iO), (26) 

where h, = (hi, .., hj, .., /i n ) and ft, 2 = hh. The transform £u(k, q) generated by 
the EbR procedure is the well-known multidimensional discrete sine transform: 

A_w(h, q) = J v dtdx^h,g( x ) i)ii(x, t), 

&h, q (x,t) := A_( e -*(h-x+ 9 t)) = (_2i)n e -iH f[ sin(/j^). ^ 27 ^ 

1=1 

For its inversion we use: 

•n iqt n r 00 00 

vh,,(x,t) = ^-j]sin(/i^), E = / ^ E •• E ' ( 28 ) 

-^1 ' '-^n . . , J , , 

/ — 1 k,q ^ mi — 1 rrin — 1 



so that (15) yields the function i/(x, i), defined in the whole space time, which 
coincides with the solution u(x, i) of the IBV problem under scrutiny for (x, t) G 
V: 

u<**) = 7±z; E •• E n^iW/^A. ^ , (x,t)G^, 

mi=l m„=ll=l K v 7 

C7(x,t) = u(x,t). (x,t)GP. 

Equation ( p9| ) implies, due to the symmetry properties of i^h,g(x } t), that t7 
is the odd (2L)-periodic extension of the solution outside V and provides the 
discrete sine-Fourier representation of the solution: 

u(x,i) = -r^lE ft s\n(hix l )[J^ t ^i~e^A.MH + 

h,Sl=l K 

E^E n Mhm) J A<? [sinfcjfo - SjOMhj, -fcf - h, • h,-)+ (30) 

sm(k J x J )v Lj (h. J ,-tf - hj ■ hj)]}, (x,i) G V, 
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where 

oc oo oo / oo \ oc 

£=£••£> £=£•£•£ (3D 

h,S mi — 1 m n — l mi-1 \nij—l/ m n — l 

and the integral J dfcj is regularized moving the singularities kj — (irrrij / Lj), rrij G 
Z + a bit off the first Quadrant. 

The Neumann problem can be treated similarly, leading to the discrete cosine 
transform: 

n 

Vh, q (x,t) := A + (e- I ( h - x +«*)) = 2"e" i « i n cos(h t xi), 

1=1 [o2) 

iqt n OO CO 

yh, g (x,t) = 2 „ £l 1 .. £ „ ^- II cos(htxt), J2 = fd q J2 ■■ E 

Z=l k,<j K mi=— oo m„ = — oo 

and to the following discrete cosine representation of the solution: 



h,cz=i n 

2E£j E n COB(fe t X t ) / ^ ^Siy^ fftcosfcj^ - ■T j )^) j (h j , -fc| - hj ■ hj)- ( 33 ) 



3=1 h^Cljtj 

cos(kjXj)wLj(hj, -kj - hj ■ hj)]}, (x, t) e V, 

where 

oo oo oo / oc \ oo 

£=£••£> £=£ • £ ••£ w 

h,C mi — — oo m Tl — — oo hj,C mi — — oo \mj=- oo / m„— — oo 

and the integral J <ifcj is regularized as in ( |30l) . 

Using the convolution theorem, one immediately recovers from equations 
( |30| ) and (|3^) the Green's integral representation ( [TT|) of the solution, corre- 
sponding respectively to the following retarded Dirichelet and Neumann Green's 
functions: 

oo oo n 

G RD (x,t ] x',t') = 1 ^H(t-t>) E ■• E e- tt <*-*) UMhxi)Mh*tl 

rn i—l Z=l 

9 n 

G flJV (x, t; x', f) = j^H(t -if) £ ' ' E e"** ) J] cos^z;) cos(/i^). 

mi£2 m^G-2 / — 1 

In the case of semicompact domains the proper spectral transform generated 
by the EbR approach has, in general, a continuous support. For instance, in 
the Dirichelet problem for the Schrodinger equation on the n-Quadrant 

V = {*EK n : Xj > 0, j = 1, .., n}, (36) 



(35) 
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the elimination operator E = A_, which restricts the definition domain to A' = 
lZ n ®I q , leads to the continuous sine Fourier transform: 

^, 9 (x,i) = (-2i) n e-^l[sin(km), y k ,,(x,f) = +1 , = J d 9 J ^ 

f=l k,g K 7j» 

(37) 

and to the continuous multidimensional sine-Fourier representation of the solu- 
tion: 



1l(nt) — -f j(k-x+gt) A_/(k,g) r dk t fk-x-fc 2 t)A f.-fk^j. 



t J ^f^e'^^W^tk,-^, (x,t)e2>. (38) 



j=l nn -l 

From the above illustrative examples it appears that the EbR procedure in 
Fourier space is very effective and, perhaps, simpler than the Green's function 
approach, which is its counterpart in configuration space. The comparison be- 
tween these two methods of elimination in examples in which the GF approach 
fails is postponed to a subsequent paper. 

3 The Elimination-by-Restriction Approach: 
The Nonlinear Case 

In this section we turn our attention to IBV problems associated with nonlin- 
ear evolution PDE's which are integrable by the inverse scattering (spectral) 
transform method. Because of the limited scope of this paper, we content our- 
selves with illustrating our method of solution by considering, as a prototype, 
the nonlinear Schrodingcr (NLS) equation 

iqt + Qxx + c\q\ 2 q = 0, q = q(x,t), (39) 

where c is an arbitrary real parameter, but this method applies as well to other 
1 + 1 dimensional soliton equations (f.i. to the Korteweg de Vries (KdV) equa- 
tion). Moreover, we confine our treatment below to solutions of (|3^) in the first 
quadrant of the (x, t) plane, namely on the semifine < x < oo for t > 0. 
Besides the initial value q(x, 0) = qo(x), the boundary value which uniquely 
specifies the solution is 

f(t) = ai v(t) + a 2 w(t), t>0, (40) 

where we have set 

v(t)=q(0,t), w(t)=q x (0,t). (41) 

Here a\ and a 2 are given real constants and, if a 2 = (a± = 0) this is the 
Dirichclct (Neumann) BV problem. Thus the problem is that of constructing 
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the solution q(x, t) of ( |3S| ) when the initial value qo(x) and the boundary value 
f(t) are given functions in an appropriate functional space (we may assume that 
they are complex valued functions which rapidly decay as x — > oo and t —> oo). 

The key-property of the NLS equation fl39| ) is that it is the integrability 
condition for the following pair of 2 x 2 matrix linear Ordinary Differential 
Equations (ODEs), 

V x = (ika 3 + Q)S>, ^ t = 2ik 2 [a 3 ,^]+M^ (42) 

where a 3 = diag (1,-1) and 

Q(x,t)= ( ^° f) \ M(x,t,k) = 2kQ-ia 3 Q x +iQ 2 a 3 . (43) 

The solution ^(x, t, k) of these equations is defined by the asymptotic condition 

V(x, i, fc)e~ ifcE,T3 -» /, x -> 00 (44) 

which uniquely defines the scattering matrix S(k, t) in the standard way, namely 
as the boundary value 

S(k,t) = V(0,t,k). (45) 

Well-known facts, which will be instrumental in the method below, are the 
following. The matrix solution ^ and, therefore (see ([l5|)), the scattering matrix 
S have unit determinant, 

det t, k) = det S(k, t) = 1. (46) 

Moreover, the property 

Qt = -CQCT 1 , C := f J ° ) (47) 

of the matrix Q, see (^3|), induces the corresponding property 

&(x,t ) k)=CV- 1 (x ) t,k)C- 1 i S\k,t) = CS- 1 (k,t)C-\ (48) 

on the Jost solution \& and on the scattering matrix S (the superscript t indicates 
hermitian conjugation). As a consequence, it is convenient to parametrize the 
matrix S by introducing the two functions a(k,t) and /3(k,t) according to the 
definition 

w-Uij <49) 

As for the fc-dependence (here the complex spectral variable k plays the same 
role as the Fourier variable in the linear case), the functions a(k,t) and (3{k,t) 
turn out to be analytic in the UHP (im k > 0) and to have there the asymptotic 
behaviour 

a(k,t) = l + Oik- 1 ), p{k,t) ^Oik- 1 ), (50) 
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for large |fc|. Finally, we remind the reader that these analyticity properties of 
a(k,t) and f3(k,t) provide the way to solve the inverse problem, i.e. S(k,t) — > 
Q(x,t), for any fixed t > 0; the basic equations of the inverse problem, which 
are not reported here, read as either Cauchy-type integral equations in the k 
variable or, equivalently, as Marchenko-type integral equation in the x- variable. 

Let us now look at the time evolution. Here the real crux of the spectral 
method appears in the evolution equation of the scattering matrix, see ( [42| ) and 

S t = 2ik 2 [a 3) S] + Z(k,t)S, (51) 

since the matrix Z(k, t) has a separate dependence on both the boundary data 
v(t) and wit) (see (fill)) according to the following expressions 

Z(k, t) = 2kV(t) - ia 3 W(t) + iV 2 (t)a 3 , r __v 
V(t) = Q(0,t), W{t)=Q x (0,t). [0Z> 



As a consequence, the evolution equation ( J51| ) cannot be immediately integrated 
to yield the scattering matrix whose knowledge is essential to reconstruct Q(x, t) 
via the solution of the inverse problem. Since only the boundary datum (^0|) is 
given, in analogy with the elimination strategy of the linear case, we introduce 
at this point the novel matrix 

S(k,t) = A- 1 (k)S- 1 (-k,t)A(k)S(k,t), A(k) := aj + 2ika 2 <j 3 , (53) 

because of its two important properties. First, its determinant is unit (see (^)) 
and its asymptotic value as \k\ — > oo is the unit matrix: 

det S(k,t) = 1, S(k,t) = I + 0(k- 1 ); (54) 

second, it satisfies an evolution equation which contains only the given boundary 
value (|40|), namely 



S t = 2ik 2 [a 3 ,S]+4kA- 1 (k)S~ 1 {-k 1 t)F(t)S(k,t), 

-c/(t) ^ (55) 
/(*) o 



F(t) = a{V{t) + a 2 W{t) = 



Though this is an important step, it does not yield the solution of our problem 
since the unknown scattering matrix S(k, t) still appears in the evolution (|55|a). 
Thus one has to find the way to relate S(k,t) and S(k,t) to each other. The 
relation S(k, t) — > S(k, t) is of course trivial as it is given by the definition ( |53"|) 
itself. This relation yields the initial value S(k, 0) for the integration of the evo- 
lution equation (f|), i.e. Q(x_, 0) -> *(x, 0, k) -> S{k, 0) = *(0, 0, k) -> §(k, 0). 
As for the inverse relation, S(k,t) —> S(k,t), one has instead to set up a RH 
problem, which finally leads to a Cauchy-type integral equation. Starting with 
rewriting ( |53| ) in the form A(k)S(k, t) = S(—k,t)A(k)S(k,t), and noting that 
the first column of S(k, t) in J49| ) is analytic in the UHP with the asymptotic be- 
haviour (pQ), one can write down two coupled integral equations for a(k,t) and 
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(3(k, t) in terms of S(k, t) by going through the standard RH problem technique. 
By assuming, just for the sake of simplicity, that no poles occur in the UHP, 
and by substituting S(k, t) with its expression obtained by formally integrating 
the evolution equation ( |55| ) , one finally ends up with the two coupled nonlinear 
integral equations 

~°° , , 2 ( 56 ) 

— oo 

whose nonlinearity is due to the fact that the function h(k,t) depends itself on 
the unknowns a{k, t) and (3(k,t) through the integral (in t) relation 

h(k,t) — {a(fc)ao(fc)/3o(— k) — a(— fc)ao(— k)/3o(k)— 

4k J dt'e ilkH ' [f(t')a(k, t')a(-k, t') + cf(t')(3(k, t')(3(-k, If)]}/ 
o 

{a(— k)a (-k)a {k) + ca(k)[3 (—k)l3 (k)- 

t 

-Akc J dt'[f(t')a(-k,t')/3(k,t') - f(t')a(k,t')/3(-k,t')]}, 
o 



(57) 



where a(k) = a\+2ika2 and ao and /3o are, respectively, the known initial values 
a(k,0) and 0(k,O). 

The discussion of these, admittedly complicate, equations, and of their im- 
plications in various directions is beyond the scope of this paper, and will be 
reported elsewhere. Here we merely note that this formulation naturally sin- 
gle out the so-called linearizable IBV problems, these being those for which 
the boundary value f(t), see (f!o|), vanishes: f(t) = 0. Indeed, in this case, 
the kernel function h(k,t) (pq ) does not depend on the unknowns a(k,t) and 
(3(k,t), and the equations (pm become linear. As a side remark, we also note 
that setting c = eliminates the nonlinearity in all the formulae given above, 
so that ([39]) becomes the linear Schrodinger equation and our equations (p6|), 
( p7| ) yield a(k,t) — 1 while (3(k,t) coincides with the usual explicit expression 
of the Fourier transform of the solution q(x,t). 

Finally, we deem of interest to report here also an approach to the IBV 
problem for the NLS equation which is different from the one given above, 



and yet equivalent as it eventually leads to the same equations (56) and fl57|). 
The main feature of this approach is that the IBV problem is reformulated on 
the whole line, i.e. for x £ (—00, 00), and that the matrix S(k, t), as defined by 
(p3[), acquires now a spectral meaning within the standard direct and inverse 
problem associated with the Lax equation, say the first ODE in (|4^). The price 
one pays to arrive at this more familiar formulation is that the nonlinear PDE 
one has to solve now is the NLS equation with a inhomogeneous source term 
rather than the NLS equation (B9) . This approach is briefly sketched here with 
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two limitations which are merely dictated by the sake of simplicity; namely we 
confine our treatment to the Dirichelet and Neumann BV problems and, second, 
we assume that the spectral data at any time t > have no discrete spectrum 
component. 

The starting observation is that, if q(x, t) is a solution of the NLS equation 
( |39|) for x € (0, 00) and t > 0, then the function 

q(x,t) =q(x,t)H(x)-r)q(-x,t)H(-x), 77 = ±1, (58) 

as defined for any real value of x, satisfies the PDE 

iqt + qxx + 2c\q\ 2 q = (1 + ^(^'(a:) + (1 - T])w(t)5(x). (59) 

where 5(x) is the Dirac delta distribution, 5'(x) is its derivative and v, w are 
defined in (fill). Obviously, 77 = 1 (77 = — 1) is the appropriate choice when one 
deals with the Dirichelet (Neumann) IBV problem. 

As implied by the spectral method based on the Lax equations, it is conve- 
nient to rewrite ([58|) and (|5^) in matrix form by introducing the 2x2 off-diagonal 
matrix (see (f43|)) 

Q(x,t)=Q(x,t)H(x)-r)Q(-x,t)H(-x), 77 = ±1, (60) 

and the PDE 

iQt-<r 3 (Q xx ~2Q 3 ) = i:(x,t) (61) 
which is, of course, equivalent to ( |5^ ) if the source term is (sec (p2|b)) 

£(z, t) = -a 3 [{l + v)V(t)6'(x) + (1 - r,)W(t)5(x)]. (62) 

The spectral approach to the equation ( |6l| ) is based on the spectral equation 

V x = (ika s + Q(x,t))V, 9 = V(x,t,k), (63) 

and it is standard. The Jost solution ^ is defined by the asymptotic condition 
([i"4"|), ^/exp(— ikxa^) — > I, x — > 00, which readly provides its expression in terms 
of the solution ^(x,t, k) introduced above on the semiline, 

*(x, t, k) = ^{x, t, k)H(x) + EV(-x, t, -k)ES{k, t)H(-x), (64) 

where E = diag (1,77) and 

S(k,t) = ES- 1 (-k,t)ES(k,t), (65) 

is precisely the scattering matrix which is defined in the usual way, namely 

^(x,t,k) -> e lkxa3 S(k,t), x^-00. (66) 
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At this point we note that this scattering matrix S(k, t) coincides with the 
matrix (|5^) with an = in the Dirichelet case (r\ — 1) and with a\ = in the 
Neumann case (77 = —1). 

It is common expedient now to introduce also the other Jost solution of ( p3| ) , 

$(x,t,k) = $(x,t,k)S- 1 (k,t), (67) 

and to take into account the identity 

00 

S t + 2ik 2 {S,a 3 } =1 [ dx$- 1 (x,t,k)[iQ t -a 3 (Q xx -2Q 3 )}q(x 1 t,k), (68) 



which, together with the inhomogeneous PDE (|6l|), entails the evolution equa- 
tion for the scattering matrix, 

00 

S t = 2ik 2 [S,a 3 ] +i J da;$- 1 (x,t,fc)I](a;,t)*(a;,i,A:). (69) 
—00 

It is now easy to show that inserting in the integral in the RHS of this equation 
the expressions (|67|), (f3J) and ( |62|) yields precisely the evolution equation ( jiiBT) 
for the Dirichelet and Neumann IBV problems, say 

S t = 2ik 2 [a 3 , S}+2k(l+ v )S- 1 {-k, t)V(t)S(k, t)-i{l-ri)a 3 S- l {-k, t)W(t)S(k, t). 

(70) 

We end this paper remarking that a good side of the present approach, which 
we will refer to as the "source-method" , is that one may take advantage of the 
more traditional inverse scattering (spectral) technique on the whole line. In 
particular, to investigate the large time behaviour of the solution q(x,t) of the 
IBV problem since the asymptotic expression, if the boundary value rapidly 
vanishes as t —> 00, are readly at hand in the usual spectral theory on the whole 
line. 



4 Literature 

The classical idea of eliminating unknown boundary values from the represen- 
tation of the solution is the essence of the Green's function approach jO, which 
makes essential use of the image method || to construct the proper Green's 
function which eliminates the unknown boundary values. The Elimination-by- 
Restriction approach introduced here is a natural and effective implementation 
of the elimination strategy in Fourier space. To the best of our knowledge this 
method for linear PDEs has never been presented before. 

An alternative method, which we term the Analyticity approach, is also 
possible, and it is motivated by Fokas discovery of the global relation and of 
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its use to solve IBV problems ||, 0], [0], 0. Our contribution to this 
method consists in using systematically the analyticity properties of all the 
Fourier transforms involved in (||), to derive a set of analyticity constraints 
which allow one to express unknown boundary values in terms of known ones 
and, in general, to study the unique solvability of IBV problems (see @]). 

Different approaches to deal with the problem of unknown boundary data 
in the study of IBV problems for soliton equations have been developed during 
the last few decades. In ||, Fokas introduced a nonlinear analogue of the sine 
transform. In j[o|, Sabatier constructed an "elbow scattering" in the (x,t)- 
plane to deal with the semiline problem for KdV, leading to a Gel'fand - Levitan 
- Marchenko formulation. In [|[ || a different approach, based on a simultaneous 
x-t spectral transform, has been introduced by Fokas and rigorously developed 
in 1 11, 1^|, to solve IBV problems for soliton equations on the semiline. It allows 
one for a rigorous asymptotics [fl3| and captures in a natural way the known 
cases [fl4| of linearizable boundary value problems. In [ fl5| we have introduced 
two alternative approaches to the study of IBV problems for soliton equations 
on the segment and on the semiline. In the first method we expressed the 
unknown boundary values in terms of elements of the scattering matrix S(k, t), 
thus obtaining a nonlinear integro-differential evolution equation for S. In the 
second method, which can be viewed as the nonlinear analogue of the EbR 
approach developed in §2, we constructed the nonlinear evolution equation ( |55| ) 
for S, which does not contain unknown boundary values and captures in a 
natural way the case of linearizable IBV problems. 

In some nongeneric cases of soliton equations corresponding to singular dis- 
persion relations, like the stimulated Raman scattering (SRS) equations and the 
sine Gordon (SG) equation in light-cone coordinates, the evolution equation of 
the scattering matrix does not contain unknown boundary data. The SG equa- 
tion on the semiline has been treated using the x-t spectral transform p); the 
SRS and the SG equations on the semiline have also been treated using a more 
traditional x- transform method respectively in Jl^l by Leon and Mikhailov and 
in [jl7| by Leon and Spire; the x- spectral data used in this last approach satisfy 
a nonlinear evolution equation of Riccati type. 

Apart from the simultaneous x — t transform, all the above approaches are 
based on the traditional 1ST [j^]. The spectral formalism for studying forced 
soliton equations has been developed by several authors, expecially in connection 
to the theory of perturbations (see, for instance, |l9[ ). 

The results presented here in the nonlinear context are: i) the construction, 
via the RH technique, of the closed form integral equation (p^), ([57]) satisfied 
by the elements of the scattering matrix, ii) The formulation of the Dirichelet 
and Neumann IBV problems on the semiline for soliton equations as forced 
initial value problems on the whole line. The equivalence between the semiline 
and whole- line problems has been already used in M, although the relevant 
equations for the spectral data given there differ from those presented here. 
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